PERIODIC BINARY HARMONIC FUNCTIONS ON LATTICES 



MIKHAIL ZAIDENBERG 

Abstract. A function on a (generally infinite) graph F with values in a field K of 
characteristic 2 will be called harmonic if its value at every vertex of T is the sum of its 
values over all adjacent vertices. We consider binary pluri-periodic harmonic functions 
/ : F2 = GF(2) on integer lattices, and address the problem of describing the 

set of possible multi-periods n = (rii, . . . ,ns) € N** of such functions. Actually this 
problem arises in the theory of cellular automata jMOWl ISulL ISu4[ IGKW| . It occurs 
to be equivalent to determining, for a certain affine algebraic hypersurface Va in 
Ap^, the torsion multi-orders of the points on Vs in the multiplicative group (Fj 
In particular V2 is an elliptic cubic curve. In this special case we provide a more 
thorough treatment. A major part of the paper is devoted to a survey of the subject. 

In mathematics, our role is more that of servant than master. 
Charles Hermite, from reminiscences by Jacques Hadamard. 
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Introduction 

We use the following notation. We let T = (vert(r), edg(r)) be a (not necessarily 
finite) graph without loops and multiple edges, i^' be a field of characteristic 2, k = 
¥2 = GF(2) the binary field; moreover, we write k = F2 for its algebraic closure. 

For i^'- valued functions / : vert(r) K, we consider two kinds of laplacians: Ap is 
averaging over balls of radius 1, respectively Ap is averaging over spheres of radius 1 
(cf. e.g., jCat IBap so that 

A+ = id+Ar, where A^{f){v) = ^ /(^/) , t; e vert(r) . 

[t),t)']e edg(r) 

Actually A^ is the a^- cellular automaton studied e.g. by Martin, Odlyzko and Wolfram 
|MUWj ■ Sutner |Sulj - [5u4] . Goldwasser, Klostermeyer and Ware |GKWj . Sarkar, Barua 
and Ramakrishnan |BR| ISBj . Hunzikel, Machiavello and Park |HMPj e.a. 

A function / on vert(r) is called harmonic {antiharmonic) if Ajt (/) = (Af (/) = 0, 
respectively). Thus for / harmonic, f{v) is equal to the sum of values of / over the 
neighbors of v, whereas for / antiharmonic this sum is always zero. We let 

HarmJ(r) = ker (Ajt) and Harm^(r) = ker (Af ) = ker (id+Ajt) 

be the corresponding subspaces of the vector space 3^{T, K) of all i^'-valued functions on 
vert(r). We simply write Harm^(r) = Harm^(r) when dealing with binary functions. 

The support = N^{f) = supp(/) of a nonzero binary harmonic (antiharmonic) 
function / will be called a nucleus {antinucleus, respectively) of F. Note that the 
constant function 1 is harmonic on an odd graph and antiharmonic on an even one. 
More generally, the (anti)nuclei can be characterized by the following two properties: 

• Every nucleus A^"*" of F is an odd subgraph, that is each vertex of A^"*" is of odd 
degree within A^"*". Whereas every antinucleus is an even subgraph. 

• Every vertex v G vert(F) \ vert(A^) has an even number of neighbors in A^^. 

Nuclei are even-parity subgraphs of F, see e.g. Amin, Slater and Zhang |ASZj . Kloster- 
meyer jKj and literature therein on parity domination in graphs. The set of all 
(anti)nuclei in F is closed under symmetric difference. 

Definition 0.1. We say that F is harmonic {antiharmonic, respectively) if there exists 
a nonzero binary harmonic (antiharmonic, respectively) function on F. A biharmonic 
graph is a graph that is simultaneously harmonic and antiharmonic. 

A locally finite graph F will be called even {odd, respectively) if the degree of every 
vertex f of F is even (odd, respectively). Every odd graph is harmonic and every even 
one is antiharmonic. Moreover, any (anti)harmonic graph can be obtained from an 
odd (respectively, an even) one by adjoining a certain number of new vertices, each 
one being joined with the old ones by an even number of edges. 

For instance, an infinite plane hexagonal lattice is odd and therefore harmonic, 
whereas an infinite plane triangular lattice is even and so antiharmonic. Actually 
the latter one is biharmonic, the nuclei being the maximal inscribed hexagonal lat- 
tices. Similarly every integer lattice is biharmonic when regarded as a graph with 
edges parallel to the coordinate axes. In all these examples, none of the (anti) nuclei is 
finite. However it is easy to find infinite graphs with finite nuclei. 
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We let spec^(r) be the spectrum of the laplacian Ap in the algebraic closure k. 
Since the laplacian Ap is defined over the binary field k, we have 

:= dimfc(Harm^(r)) = dim^(Harm^(r)) . 

Thus by definition, F is harmonic if and only if G spec+(r), antiharmonic if and only 
if G spec~(r), and biharmonic if and only if 0, 1 G spec^(r). The shift by 1 being an 
involution on k, this justifies our terminology. 

Our aim is to determine the set of all multi- indices n = (ni, . . . , Ug) G N'^ such that 
the integer lattice possesses an n-periodic nonzero binary harmonic function. For 
instance, for n = (1, . . . , 1) such a function / on Z'' must be constant = 1, but the 
constant function 1 on Z'^ is not harmonic, although it is antiharmonic. 

Given a Galois covering vr : F' ^ F with the Galois group G, there is an isomorphism 

TT* : Harm^(F) ^ [Harm^(F')]^ , / i — ^ / o tt, 

where the latter vector space consists of all G-stable (anti) harmonic functions on F'. 
Thus our problem reduces to the following one: 

Letting C„ be a circular graph with n vertices, we consider a finite abelian group 
Ijn := Y[i=i '^ln{L and the associated toric lattice (the Caley graph of Z^j) 

Tf^ = C„^ X . . . X C„^, n= (ni, . . . , n,) G . 

Then J/ can be viewed as the maximal abelian cover of with the Galois group G = 
'Y^l^^riiLei C Z*. Hence the space [Harm'''(Z^)]'^ of pluri-periodic binary harmonic 
functions on Z* with periods niCi, . . . ^UgSs can be identified with Harm"'"(T^). So we 
would like to determine the set of all harmonic toric lattices T^. 

In section 1 we deal with (anti)harmonic functions on trees. Following Amin, Slater 
and Zhang jASZj . Gravier, Mhalla and Tanner [(ilMTj we provide in 11.31 an algorithm 
that computs the dimension of the vector space Harm^(F). 

In section 2 we give an account for some spectral properties of laplacians on mul- 
tidimensional grids and toric lattices. We also mention some uniqueness theorems for 
binary harmonic functions on graphs and an application to the game 'Lights Out'. One 
of the typical results of this section is as follows fsee I2.22l b). 

Theorem 0.2. Given n = (rii, . . . ,n<j) G N'^, there exists a nonzero n-periodic binary 
harmonic function on If (that is, the toric lattice is harmonic) if and only if the 
affine variety in A| with equations 

s 

(1) ^(Xi + xri) = 1, <' = 1, t = l,...,s, 

i=l 

is nonempty. 

Section 3 is devoted to 2-dimensional grids and lattices. According to the above 
theorem and taking into account the covering trick, to distinguish the harmonic tori 
it is enough to determine all bi-torsions (ord x, ord ?/) G N^^j ^ of points {x,y) on the 
elliptic cubic curve E* in A| with equation 

X + 1/x + y + 1/y = 1. 
""^We denote by ord^ € Nodd the multiplicative order of an element ^ G A:^ . 
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As suggested by Zagier, we consider an equivalence relation on the set Nodd of all odd 
natural numbers defined by the connected components of the partnership graph 7^^^ 
with the set of vertices Nodd and the set of edges 



We indicate some simple properties of this graph. In particular, all its connected 
components are finite (Theorem IH.llj) . In Appendix 1 the reader will find an overview 
of the first 13 connected components of T^^-* calculated by Zagier with PARI. 

Finally in Appendix 2 we provide a survey on binary Chebyshev-Dickson and Fi- 
bonacci polynomials, as these are closely related to our subject. 

We are grateful to Don Zagier for his clarifying suggestions in section 3 and the 
calculations in Appendix 1. Our thanks also to Roland Bacher, Silvain Gravier, Lenny 
Makar-Limanov and Pieter Moree for highly useful discussions and pointing out refer- 
ences, to Andrzej Schinzel for disproving a conjectural inequality for the Euler function 
in I3.141 2. to Gottfried Barthel and Roland Bacher for helpful editorial remarks, and to 
Audrey Inshakov for his assistance with MAPLE. Gottfried Barthel also helped with 
some computations for the Euler function. 

The aim of the present survey being rather pedagogical, we have to apologize that the 
list of references is by no means complete, neither our survey follows the chronological 
order of events. 



The following proposition is well known, see e.g., |Anj . |Su2| 4.1-4.3], cf. also |DGj . 
|OZj . For the sake of completeness, we give a simple argument that applies in the 
characteristic 2 case. 

For a vertex v G vert(r) (respectively, for an edge [u,v] G edg(r)) we let F — f 
(respectively, F — [m,!"]) be the graph obtained from F by deleting v and all incident 
edges (respectively, the edge [u,v] but not the vertices u and v). We let adj(F) be the 
adjacency matrix of a finite graph F. Notice that adj(F) coincides with the matrix of 
the laplacian Ap. We let Xr(^) the characteristic polynomial of adj(F), and denote 
by (l)j, the number of i-matchings in F i.e., of all possible choices of i non-incident 
edges among the e edges of F. 

Proposition 1.1. For any finite graph F with n vertices and e edges, the following 
hold: 

(a) Xr(^) ~ Ylt^o^ (f)r"^"~^*' particular n and Xr (^^^ of the same parity. 



{ [ord X, ord ?/] : {x,y) E E*} . 



1. Harmonic forests 



(b) Wv G vert(F) 






[v,v']£edg(r) 



(c) y[u,v]eedg{T), 

XA^) = Xr~[u,v](^) + Xr-{u,v}ix) ■ 

(d) V[u, f ] G edg(F) with degw = 1, 

Xr {x)=x- Xr_^ {x) + Xt-{u,v} i^) and so, Xr (0) = Xr-{u,v} (0) • 
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(e) Given u,v,w G vert(r) such that degu = degv = 1 and [u,w], [v,w] G edg(r) 
(that is u,v are extremal vertices ofT joint with w) one has 

Xri^) = x^-Xr-{uM(^) ^^'^ ^r(l) = Xr-Kt,}(l) • 

(f) Given u,v,w G vert(r) such that degu = 1, degf = 2 and [u,v], [v,w] G edg(r) 
(that is [u,v] is an extremal edge ofT joint with w) one has 

Xr{x) = (1 + x^)-Xr-{u,v}^^^ + ^'Xr-i™} W (^^d so, Xr(l) = Xr-{„}(1) • 

Proof. The order n symmetric determinant det = XIo-gs ^<^' '^here as usual Sn stands 
for the n-th symmetric group and rria- = ± nr=i '^i,o-(i)' I'educes modulo 2 to ^o-2=id"^o- 
by cancelling equal terms rua and m^^i with cr^ 7^ id. This leads to (a). Now (b) and 
(c) can be easily deduced from (a). In turn (d), (e) and (f) can be deduced by virtue 
of (b). □ 

1.2. Let r be a finite forest that is, a disjoint union of trees. It can be reduced, in two 
different ways, to a rather simple one, by 

• iteratively suppressing an extremal vertex (leaf) as in (d). In this way we finally 
reduce F to a forest r~j with only isolated vertices; 

• iteratively suppressing a pair of extremal vertices as in (e) or a pair of extremal 
edges as in (f). Via this procedure, T will be finally reduced to a forest 
with only isolated vertices and isolated edges. 

This gives the following result, see p\SZj or, in any positive characteristic, [G MTt 
Theorem 4 and Corollary 6]. We recall the notation d^(r) = dim(Harm^(r)). 

Corollary 1.3. (a) A forest T is harmonic (antiharmonic, respectively) if and only 
'^f^red O^red' respcctivcly) contains an isolated edge (an isolated vertex, respec- 
tively). 

(b) Moreover, for any (^red' respectively), the number of isolated edges (of 
isolated vertices, respectively) is d'^(T) (d"(T), respectively). 

Proof, (a) For a disjoint union F = F' U F" of two graphs we have Xr — Xr'Xr"- Since 
det(A];) = xr(0) and det(A^) = Xr(l), we get det(A^) = det(A^,) det(Ar„)- By 
virtue of (d)-(f), the first reduction preserves Xr(0)) second one Xr(l); that 

det(Ap) = det(A^± ). Thus det(Ap) = 1 if and only if F~^ is empty, and det(Ap) = 1 

red 

if and only if F^^ consists of isolated vertices. This proves (a). 

(b) Following our iterative procedure, we can easily see that every (anti) harmonic 
function on F restricts to a (anti)harmonic function on F^^ (on F~^, respectively). 
Moreover we can reconstruct the (anti)harmonic functions on F from their restrictions 
to F^^, respectively. Indeed, for any isolated vertex v of F~^, the 5-function 5y on 
vert(F~j), which takes value 1 at w and at any other vertex, is antiharmonic. At 
every step, Sy uniquely extends from a smaller graph to a bigger one preserving anti- 
harmonicity. This results finally in an antiharmonic function Sy onT. 

On the other hand, given an antiharmonic function /i on F, it is uniquely determined 
by the restriction h \ F~j. This restriction can be decomposed in the basis of 6- 

functions (5^ : f G vert(F~^)) in 5'(F~j,fc). Hence (^6y : v G vert(F~^)j form a basis 

of Harm~(F). 
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Similarly, given an isolated edge [u^v] of F^^, 5[u,v] = 5« + 5^, is a harmonic function 
on r^^. At each step it extends uniquely to a harmonic function on a bigger graph, 

G Harm'''(r). These functions form a basis of Harm^(r). 

□ 



and finally to a function 5^.^^., 
This shows (b). 



Remark 1.4. The analysis of (anti)harmonicity of unicyclic graphs can be reduced in 
the same way to that of cyclic graphs C„ |Su2l §4]. As for the latter one, see section 
2.1 below. 

2. Chebyshev-Dickson-Fibonacci polynomials and harmonicity 

2.1. 1-dimensional case. We refer the reader to Appendix 2 for a survey on the 
Chebyshev-Dickson polynomials T„ {En) of the first (second) kind and the Fibonacci 
polynomials F„. We also need the following notation. 

2.1. For n G Nodd, the order and the suborder of 2 modulo n are, respectively, 

f{n) = ord„2 = min{j : 2^ = 1 mod n} 

and 

fo{n) = sord„2 = min{j : 2^ = ±1 mod n} . 

Thus /(n)//o(n) G {1,2}. Moreover, 

f(n) = 2fo{n) 3j G N : 2^ = -1 mod n . 

Letting q = 2^°^"'^ we note that n \ (g — 1) if fo{n) = f{n) and n \ (g + 1) otherwise. 
Anyhow, n divides exactly one of g — 1 and g + 1. Further, /(2^' — 1) = /o(2^' — 1) = r 
Vr > 3 (but /o(3) = 1, /(3) = 2) and /o(2'' + l) = r = /(2^ + l)/2 Vr > 1, see Appendix 
Bin IMDW] . 

2.2. We notice that A^^ = r + r~^, where r G End(A^) is the cyclic right shift, and 
Ap^ = Ti + Tr, where ti (r^) G End(A^) is the left (right) shift. Hence the adjacency 
matrices of the graphs C„ and P„ are, respectively, 

/010...00l\ /010...000\ 

1 1 ... 



adj(C„) 




1 



1 


























1 







1 

0/ 



adj(P.) 



1 



■•. 10 
... 1 1 
\ ... 1 / 

En, respectively. For n odd 







with the characteristic polynomials Xc„ — ^« Xp„ 
we have spec(r) = fin, where /i„ stands for the cyclic group of n-th roots of unity in 
k. According to the spectral mapping theorem, spec(A^^) = + ^"^ : ^ G fin} and 
spec(Aj ) = {1 + ^ + '■ ^ £ f^n}- Moreover, for the circular graphs C„ the following 
results hold, see e.g. [Eu3l 4.1, 6.1], jEHl 2.1], [Hj, 3.3.8]. 

Proposition 2.3. (a) Vn > 3, C„ is antiharmonic (that is A^^ is non-invertible). 
Whereas Cn is harmonic (i.e., is non-invertible) if and only ifn = Omod 3. 
(b) The minimal polynomial of is Fk if n = 2k and xRk, where Rk = VF^, if 
n = 2k + l. 



PERIODIC BINARY HARMONIC FUNCTIONS 7 

(c) The polynomial xRk having simple roots, for every n G Nodd ^he matrix adj (C„) 
is similar over k to the diagonal matrix diag(C* + : i = 0, . . . ,n) , where 
( & fin is a primitive n-th root of unity. 

(d) Consequently, Vn G Nodd, {^0^ = ^c„' ^^^^^ 1 = 

(e) The kernel of A'^_^^ is two-dimensional, spanned by the vector (1, 1, 0, 1, 1,0,.. .) 
and its shift. If n is even then the kernel of is also two-dimensional, 
spanned by the vector (1, 0, 1, 0, 1,0,.. .) and its shift. For n odd this kernel is 
one- dimensional, spanned by (1,1,1,1,1,...). Hence d'^{C^k) = 2 = d~{C2k) 
and d-{C2k-i) = 1 VA; > 1. 

(f) Respectively, the nuclei ofC^k are the cyclic shifts of = {vi : i ^ mod 3}, 
the antinuclei of C2k are the cyclic shifts of N" = {vi : z ^ mod 2], whereas 
N" = Cn is the only antinucleus of Cn, Vn = 2A; — 1. 

Similarly, for the paths P„ we have the following results, see e.g. |M()Wj . |Su3j . |BR| 
4.4], |Ml 3.3-3.4]. 

Proposition 2.4. (a) 3(Ap^_J"^ n G Nodd, and 3(A^^_J"^ 

n^O mod 3. 

(b) The minimal polynomial of Ap^ is En- 

(c) Vn G Nodd, Ap^ admits a generalized inverse np^ G End(5'(P„, /c)) such that 
Ap Kp^Ap = Ap . 

(d) Vn G Nodd, oi'd(Ap^_J = 2e„ - 2, where e„ = min{j G N : (A^^)^' = A^^}, 
Cn e Neven is such that (e„ - 1) I (g - 1) for q = 2^o("). 

(e) Vra G Nodd, (AJ^ J^"? = (A^^ J^. Furthermore, if n e Nodd and n ^ Omod 3 
then ord(Aj^_J | (2g - 2), where q = 2^o(n). 

(f) The only nucleus of P^k-i is = {vi : i ^ mod 3}, and the only antinu- 
cleus of P21-1 is = {vi : i = 1 mod 2}. Hence d^{P3k-i) = d^{P2i-i) = 1 
VA;,/ G A^. 

Remark 2.5. For every n G Nodd, both (Ap^ and the generalized inverse Kp^ for 
Ap^ are explicitly found in |SBj . 

2.2. Spectra of products. 

2.6. Letting E, E' be vector spaces over a field K and ei, . . . , (e'^, . . . , e^, respec- 
tively) be a basis of (£", respectively), we represent every X = Xijei®e'j G E®E' 
by the matrix (or pattern) X = (xij) G Matm^niK). Following \Bl{\ §5], fSB| §4], for 
any two square matrices A G MatTO^m(-f^) and S G Mat„^„(ii') we consider the Sylvester 
derivation 

5a,b &'E'^d{E E'), X\ — >AX + XB\ with the matrix C = A iS) I + I <S) B . 

The following lemma is well known |Gal V111.3] and holds for arbitrary fields. We 
provide a simple argument in the characteristic 2 case for K = k. 

Lemma 2.7. (a) In the notation of \2. 61 we have 

Xc{.x) = Resy (xAi^ + y)^ Xsiy)) o.nd spec(C) = spec(v4) + spec(5) 

(the Minkowski sum in k '^). 

^We recall that for a commutative semigroup 11 and for two subsets Ai,A2 C n their Minkowski 
sum is Ai + A2 = {Ai + A2 : Aj e A^, z = 1, 2}. 
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(h) C is invertible if and only if the characteristic polynomials Xa '^^'^ Xb ^''^^ 
prime. 

Proof. Let A = Sa + Na be the Jordan decomposition of A G End(A^), with SaNa = 
NaSa, where Sa,Na G End(A^), Sa is semi-simple and Na is nilpotent. Then A'^ = 
S\ = Sa for certain q = 2^ , r > 0. 

We fix g = 2^ so that = Sa, = Sb and = Sc- Since A ® 1 and 1 ® 5 
commute, we have 

5c = C"^ : X I — > A'^X + X{By = SaX + XSb 

i.e., S'c = S'^ ® 1 + 1 (S> 5'^. If the bases (cj), (e^) as in 12.61 are diagonalizing for Sa, Sb, 
respectively, with SA^ei) = XiCi and Sb^g'^) = fJ^je'j, then (ei®e^) is a diagonalizing basis 
for Sc with Sc{ei ® e'j) = (Aj + fj,j)ei ® e'y For any two polynomials p = YYiLii^ + ^i) 
and g = HJlil^ + f^j) have |vdWj 

Resy{p{x + y),q{y))= JJ (x + + /xj) . 

l<i<m,l<j <n 

Since = Xs^ stc, the assertions follow easily. □ 

Remark 2.8. If C« = Sc then = p(C"2) = ^(5^) Vp G k[x]. It follows that 

Pmin(C) = 0, where 

•y=X+fi, AGspcc(A), /t6spcc(B) 

2.9. The Cartesian product F = Fi x r2 of two graphs Fi, F2 is defined via 

vert(F) = vert(Fi) x vert(F2), edg(F) = [vert(Fi) x edg(F2)] U [vert(F2) x edg(Fi)] . 

In particular, the m x n-grid is the product Pm,n = Pm x Pn, and the toric m x n-lattice 
is the product Tm,„ = Cm x C„. 

Fixing an ordering of the m (n) vertices of Fi (F2), we may regard any i^- valued 
function on Fi x F2 as an m x n-matrix X with entries in K. The laplacian Ap acts 
on X via 

: X I — > adj(Fi)± ■ X + X ■ adj(F2) = adj(Fi) • X + X ■ adj(F2)^ , 
where A' = A and yl+ = A + 1. 

We let spec^(F) = spec(Ap) C k and Xr = ^adj{r)- From ITTI we deduce such a 
corollary, see e.g. |Ba| Lemma 8]. 

Corollary 2.10. (a) The spectrum spec~(F) of the product F = Fi x F2 0/ two 
graphs is the Minkowski sum of the spectra spec~(Fj), i = 1,2. Moreover 

xAx) = ReSy {xr,ix + y),Xr2iy)) ■ 

Whereas 

spec"'"(F) = 1 + spec~(Fi) + spec~(F2) . 
(b) Consequently, F is antiharmonic if and only if the characteristic polynomi- 
Xv^i Xv2 ^'^^ coprime, and is harmonic if and only if the polynomials 
Xfi? Xr2 ^'^^ coprime. 
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2.3. 2-dimensional grids and tori. The following is immediate from 12.101 see jBRj . 
jSuSj, pMP] . 

Proposition 2.11. (a) The grid Pm-i,n~i is antiharmonic (respectively, harmonic) 
if and only if the Chebyshev- Dicks on polynomials E^-i and En-i (respectively, 
Em-i and E^_^) are not coprime. Furthermore Pm-i,n-i is antiharmonic if and 
only z/gcd(m, n) ^ 1. 

(b) det(A+^_^_J = Res,.(E„_i,E+_i). 

(c) Vm, n > 3, the toric lattice „ is antiharmonic and, moreover, is an even 
graph. Furthermore T^.n is harmonic if and only if the polynomials and 
are not coprime. 

Corollary 2.12. (a) Vfc, / G N, the grids P2k-i,3i-i and P3k-i,2i-i are harmonic. 
The grid Pm-i,n-i different from any one of these is harmonic if and only if the 
system 

(2) u + u-'^ + v + v-'^ = l = u"' = v'' 

admits a solution {u,v) G (A;^)^. 

(b) We have 

spec(AT^^^^) = {u + u-'^ + V + v~'^ : u e fim, v E , 
respectively, 

spec{Aj^ = {1 + u + + V + v^^ : u G jimi v G . 

Thus Tm,n is harmonic if and only if the system (0j admits a solution, if and 
only if either mn = mod 3 or Pm-i,n-i is harmonic. 

(c) \fm,n = mod 5, both the grid Pm-i,n-i and the toric lattice Tm,n are har- 
monic. 

(d) IfTm,n, respectively, Pm-i,n-i is harmonic then so is Tkm,in' respectively, Pkm~i,in~i 
Vfc,/GN. 

(e) Vg = 2", Vg' = 2^, T^.n is harmonic if and only if Tg^.g'n is, and Pm-i,n-i 
is harmonic if and only if either Pqm-i,q'n-i is, or one of the following holds: 
m = mod 2, n = mod 3 or m = mod 3, n = mod 2. 

(f) In particular 'i a, h > 0, Tq^qi and Pg-i,g'-i are not harmonic. 

Proof. Pm-i,n-i is harmonic if and only if Em-i{z) = En-i{z + 1) = for some z E k. 
These equations are satisfied by z = (respectively, 2; = 1) if and only if m = 
mod 2, n = mod 3 (respectively, m = mod 3, n = mod 2), see I5.4[ a.e. Suppose 
further that z ^ 0,1. Letting 

z = u + u^^, z + 1 = V + v~^, where u,v E k^ , 

by virtue of in 15. II and 15.41 a we obtain 

E^_i{z) = E^^iiu + u-^) = = E+_^{z) = En-i{v + v~^) 

^ T^{u + u~^) = m'" + M-™ = = Tn{v + v-^) = t;" + t;"" . 

This shows (a). The same argument proves (b). The assertions (c), (d) and (e) follow 
from (b) by virtue of I5.4l e and, in turn, imply (f). □ 
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In order to find all harmonic toric lattices it is enough, by virtue of l2.12l e. to restrict 
to Tm,n with {m,n) G N^^^. The following facts are estabhshed in |GKWt Theorem 
14], see also PMFI 5.1]. 

Proposition 2.13. Vg = 2^, r > 1, the toric lattices Tg_i.g_i, Tq__i^q_|_i and T^+i^g+i 
are harmonic except for Ti^i and T7.7. 

Proof. Letting, according to 15. 21 c. d. 

Ag = roots(T,_i) = {0} U{ze ¥^ : Ttw,{z-') = 0} , 

Bg = roots(r,+i) = {0} U{ze ¥^ : Ttw^{z-') = 1} 

and 

A+ = roots(r+ 1) = 1 + Ag, S+ = roots(r++i) = 1 + 5, 

we have Ag U Bg = ¥g and Ag n Bg = {0}. Indeed by Og, Tg+i + Tg^i = x^+^ 
and Tg^iTg_i = x^'^ + = x'^{x'^~^ + 1)^. So the zeros of the product Tg^iTg_i fill 
in Fg, while is the only common zero of T^+i and Tq_i. Hence card(y4g) = g/2 and 
card(Sg) = g/2 + 1. It follows that Ag^B^ ^ and BgHB^ ^ 0. Thus the polynomials 
Tq_i, Tq"*;,.!, respectively, Tg+i, T^^^^ are not coprime. In view of l2.111 c. the toric lattices 
and T^+i^^+i are harmonic Vg = 2^', r > 1. 
Suppose further that r > 2 and Tq_i^q_i is not harmonic, that is Ag fl = 0. Then 
Ag Bg\ {0}. Actually A^ = Bq \ {0} as these sets have the same cardinality. Thus 
roots(T,"ti) = roots(Fg_|_i). More precisely, 

= (x+l)T+_i ^ = (x+1)2f+_i ^ xHl = . 

For every 2; G Fg\F2 we obtain z+1 = Fg_i{z+ 1) + 2;^F,„i(z). Equivalently, by virtue 
ofOld. 

(3) (z + 1) (1 + Trip,((^ + 1)-^)) = ^^Trip,(;.-^) . 

From (j2)) we deduce the following alternative. 

• Either 

= z + l =^ z e¥s\¥2 G¥g\¥2 =^ r = mod 3 , 
and then 1 + Ti^^{{z + 1)-^) = Ti^^iz-^), 

• or 1 + TrF,((;z + 1)"^) = Trp^iz-^) = and so, Fg^i{z) = G F, \ Fg. 
Henceforth, if F, D Fg and F, ^ Fg then 

card(A, \ {0}) =g/2-l>g-8 =^ g<14 =^ g = 8, 

which is a contradiction. If F, 2 ^8 then, by the same argument as above, Fg_i{z) = 
0Wze¥g\¥2. Hence 

card(Ag \ {0}) = g/2-l>g-2 =^ g<2, 

which again gives a contradiction. 

Therefore ¥g = Fg. Indeed, for g = 2^ we have = Bg\ {0} and so, the toric 
lattice Tg_i^g_i = T^j is not harmonic, as stated. □ 
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Remarks 2.14. 1. By virtue of lTT^ b. 

r = mod 2 g - 1 = mod 3 G n n 5^ , 

hence both Tg_i^g_i and Tg„i^g+i are harmonic, and 

r = 1 mod 2 g + 1 = mod 3 G A, n fi^ H fi+ , 

hence both Tq_i_g+i and T^+i^^+i are harmonic. 

2. The polynomials hi{x) = + x + 1 and h2{x) = x'^+x + 1 satisfy hi{x+l) = hi{x), 
i = 1,2. They divide the Fibonacci polynomials Fg±i (and hence also Tg±i) in the 
following cases: 

hi I -Fg-i <^=^ r = mod 4 and hi \ Fg^i <^=^ r = 2 mod 4 , 
^2 I -fg-i ^^=^ r = mod 8 and h2 \ Fg^i <^=^ r = 4 mod 8 . 

2.15. The above theory can be naturally extended to the laplacians Ap on 3'(r,Z) 
and on 3^(r, ¥p) for all primes p>2, see e.g., |M()Wj . |(4MTj . [HMP] . We say that T is 
p- (anti)harmonic if ker(Ap) has a positive dimension in 5'(r, Fp). For s-dimensional 
grids, and especially for 2-dimensional square grids, the following is proved in (HMPl 
§§4-5]. 

Proposition 2.16. (a) For n = (ni, . . . , n^) G N'^, the grid Pn is p-harmonic if 
and only if 



det(A+J= n l-E^2k+i)+^2Z%i) =0 modp, 




where c:„ = 

(b) Vn > 3 there exists a prime p such that the square grid Pn-i,n~i ond the toric 
lattice Tn^n ore p-harmonic. 

(c) Pn-i,n-i C^n,n, respectively) is p-harmonic for every prime p if and only ifn = 
mod 5 or n = mod 6 (n = mod 5 or n = mod 3, respectively). 

(d) If I > 5 and p are primes such that p is a primitive root modulo I then both the 
square grid P;_i^;_i and the toric lattice T/ / are not p-harmonic. 

(e) For every prime p with at most two exceptions, the set Ip of all primes I such that 
the square grid Pi^i^i^i (the toric lattice Ti^i, respectively) is not p-harmonic, is 
infinite. 

(f) The square grid Pn-i,n-i o-nd the toric lattice T„^„ with n = {p ±l)/2 are 
p-harmonic for every prime p > 23. 

For the proof of (a), (b), (d), (f) see |HMPj . 4.6, 4.4, 4.3, 4.7 and 5.4, respectively. 
The proof of (e) in |HMPj is based on (d) and on a result of Heath-Brown jiHB] on 
Artin's conjecture of primitive roots. ^ 



''The above product being an integer. 

^This conjecture suggests that every integer n ^ —I which is not a square, is a primitive root 
modulo I for an infinite set, say, /„ of primes I. The result of Heath-Brown loc.cit. says that the 
property in Artin's conjecture holds for all primes with at most 2 exceptions, and for all square-free 
integers with at most 3 exceptions. However, so far no concrete example of a prime satisfying the 
conjecture has been found, see | HB[IMollMu| . 
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2.17. For a graph T we let as before d^{T) = dim(ker(Ap )) = dim(Harm^(r)). The 
kernels of the grid laplacians Ap admit the following description, see |Su2j . |Su8[ 
3.6, 3.10], [EHl 3.6, 4.1] andHnni'below. 

Proposition 2.18. (a) Vp G k[x], dim(ker(p(Ap ))) = deg(gcd(p, 

(b) ker(A±^ J = ker(E^(A±J) = keT{E^{A^jy 

(c) Furthermore, ker(Ap^_^^_J ^ ker(Egcd(m,„)-i(Ap^_J). 

(d) Consequently, d-{Pm-i,n-i) = gcd(m,n)-l and d+{Pm-i,n-i) = deg(gcd(E„_i, 

(e) Vm = 2^ , Wn = 2^p, where r > 1 and p G Nodd; one has ^ 

if p ^ mod 3, {and so Pm-i,n~i is not harmonic) 
rf+(P„_i,„_i) = d+(Pn-i) = { 2^'+i if p = Q mod 3 anc? A; < r - 1, 

otherwise . 

(f) Moreover min{n > m : (i"''(Pm-i,n) = — 1} = |m — 1. 

Examples 2.19. The path P2 and the grids P2,2n-i, n > 1, P2 x Cn and Ta^^, n > 3, 
are harmonic, whereas P2,2n, n> 1, are not. The grid P2,3 has the nuclei 

1 1 1\ /O 1 0\ A 1 

1 Oy) ' Vl 1 V 1 

Thus c?+(P2^3) = 2. Similarly, the grid P2,2n-i (-P2,4n-i, respectively) has a nucleus 

1 1 1 ... 1\ /O 1 1 ... 

101010... lyl' respectively, (^^ ^ ^ g 1 1 1 . . . 1 

2.4. n-dimensional case. From ITTTll we deduce by recursion the following, cf. |SBj . 

Proposition 2.20. (a) spec"(ni=i Tj) = Z]i=i spec"(ri) anc? spec+ ([1^=1 Tj) = 
1 + Ei=ispec"(r,). 

(b) For the product T = Y[i=i of si harmonic and s — Si antiharmonic graphs 
we have {si + 1) mod 2 G spec^(r). Consequently, T is harmonic if Si is odd 
and r is antiharmonic otherwise. If at least one of the factors Fj is biharmonic 
then so is T. 

(c) V/i G Harm^(Fi), i = l,...,si, Wgj G Harm^(Fj), j = Si + 1, . . . , s, the func- 
tion h = {<S)iLi fi) ® ((S)j=si+i 9j) ^ k) is harmonic for Si odd and anti- 
harmonic for si even. 

(d) // is a nucleus of V-i, i = 1, . . . , si, and N~ is an antinucleus of Tj, j 



Si + 1,. . .,s, then N = Hiii x 11^=51+1 « nucleus ofV = [] ■=! f 

Si odd and an antinucleus of F for si even. 

2.21. We keep the notation 



or 



Tn = Y\Cn,, = n ^^"^ n - 1 = (ni - 1, . . . , n, 

i=l i=l 

For the next results see e.g. [Sul], jBH §5-6], [HMPl §3]). 



^See also Pu3l 5.2] for the case m = 3 • 2'', n = 2^p. 



PERIODIC BINARY HARMONIC FUNCTIONS 



13 



Proposition 2.22. (a) For any graph T and for every n G N, 

Xrxp„_i(a;) = ReSy{xr{x + y),Fn{y)) , Xrxc„_i(a;) = 'R-eSy {xr{x + y),Tn{y)) 
and 

spec^(r X C„) = {A G A;^ : A + A"^ G spec^(r)} . 

Hence TxCn is harmonic if and only if 1 + X + G spec+(r) for some X E fin- 

(b) We have 

s 

spec+(T^) = {1 + + C') ■■ ^i^fin,, t = l,...,s}. 

1=1 

Thus Tfj is harmonic if and only if the system (QJj in \U.^ has a solution {xi, . . . ,Xs) G 

(c) IfTxCn is harmonic then so is T x Cin for every I >1. 

(d) // Tn is harmonic then so is T^i x T^, \lfh G N*, Vn' = (/ini, . . . , IgUs) G N*, 
where Zj > 1 Vz = 1, . . . , s. 

(e) IfV X C2n (n > 3) is harmonic then so isT x Cn- Consequently, V x (r >2) 
is harmonic if and only if so is T. 

(f) = (2^^, . . . , 2'"''), the toric lattice Tn is not harmonic. 

Proof, (a) follows from 12.21 and I2.1(JI and implies (b) by recursion. The covering F x 
Cin — r X C„ with the Galois group Z//Z induces the injections 

TT* : Harm^(r x C„) ^ Harm^(r x Qn) , / 1 — > f o n . 

This proves the harmonicity of Tn> in (d), whereas that of the product x follows 
from 12. 201 b. The proof of (e) uses (a) and the fact that (j) : k ^ k, x 1 — > is an 
automorphism, (f) follows from (e) by recursion. □ 

Examples 2.23. For any antiharmonic graph F and Vn > 3, V/ > 1, the products 
F X C„ and F x P21-1 are antiharmonic , whereas F x C3; and F x P31-1 are harmonic. 
If F is harmonic then so are the products F x C„ and F x P21-1, whereas F x C31 
and F x P^i-i are antiharmonic. See also |SB|, §§5,6] for the (anti)harmonicity of the 
hypercubic grids Pn and of the products Pn x T^- 

2.5. Symmetrization. 

2.24. Let K he a field with char(i^) =2. If a : F — > F is an involution then for any 
nonzero / G 9^(F, K), either f o a = f or the average g = f + foa is again nonzero 
and is a-stable: go a = g. Anyhow, if Harm^(F) 7^ {0} then also [Harm^(F)]° 7^ {0}. 
Moreover, if F = Fix(a) ^ then for any / G [3'(F, K)]", 

(4) A^(/)|F = A^(/|F), 

and so the restriction of an a-symmetric (anti)harmonic function / to F is again 
(anti)harmonic. Furthermore, if f\F = then also /|(F F) G Harm^(F F). Thus 
if F is (anti) harmonic then so is at least one of the graphs F and F F. In case that F 
is a 'separation wall' for F, from the above discussion we deduce the following result. 

Lemma 2.25. Let a be an involution of F such that F = Fix(a) separates F that is, 
F F = F+ U F^, where F+ and F~ are two disjoint subgraphs o/F with a(F^) = F^. 
IfT is (anti) harmonic then so is at least one of the graphs F and F^. 
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Corollary 2.26. (a) //Fx P„_i is (anti)harmonic then so is T x Cn- Vice versa, 
ifTxCn is (anti) harmonic then so is at least one of the graphs T and T x P„_i. 

(b) (cf. [SB, Q.l]) Consequently, if the grid = 11^=1 -^n—i is harmonic, where 
n = {rii, . . . ,ns) then so is the toric lattice Tn. In particular for every n = 
^2»-i — 1^ ... ^ 2'^'' — 1) the grid Pn is not harmonic. 

Vice versa, if is harmonic and 11^=1 '^i ^ mod 3 then is harmonic 
too. 

(c) IfTx Pn-i is harmonic then so is T x Pin-i VZ > 1. 

(d) IfTx P2n+i is (anti) harmonic then so is at least one of the graphs T and F x P„. 

Proof, (a) follows from 12.251 To show (b), letting / be a nonzero (anti) harmonic 
function on F x Pn-i symmetric w.r.t. the reflection 

« : Pn-l — ^ Pn~l, Vi I > f(n-l-i) mod (n-1) 5 

the extension of / by zero to F x C„ ^ F x P„_i is again (anti) harmonic, as required. 
The converse in (b) follows from 12.2^ Iterating this argument yields the first and the 
last assertions of (c). The second one follows bv l2.22l d. 

To show (d) we take I copies P^*ii, i = 1, . . . , / of Pn-i- For a nucleus of F x P„_i, 
we consider its copy A^i in F x P^L\, the mirror image N2 of A^^i in F x P^^\, the mirror 

image N3 of A''2 in F x P„_i, etc. Taking also new vertices v[, . . . , v'i_i and representing 
F X Pin-i as ordered 'connected sum' of the graphs 

FxpS, FxK}, FxpS, Tx{v',}, FxK_J, F x P^^ 

we obtain a nucleus A^' = IJi=i of F x P;„_i. □ 

Remarks 2.27. 1. Starting with N = P2 the proof of (d) gives a nucleus of P31-1 
V/ > 1. 

2. Instead of taking average of / over an involution, one can consider the average of / 
over the shifts on Zn. Suppose that for some i E {1, . . . , s}, the subgroup Zf^(i) of is 
not harmonic, where n*^*^ = (ni, . . . , nj_i, n^+i, . . . , Ug). Then for any / G Harm'''(Zfi), 
the average f + fg. + . . . + of / over the shifts by the subgroup Z„. C Zfj must 

be zero. This means that the intersection of each nucleus A^"*" of with every 'line' 
Ip^i = {p, p + Ci, . . . , p + {ui — l)ej}, p G Zfi, has even cardinality. 

2.28. For a vertex f of a graph F we let 

(5) a+ = 6,+ Yl Sue3^{r,k). 

[u,v]eedg{T) 

Actually the toric lattice represents the Caley graph of the group Zn = 111=1 '^/^iZ 
with its standard generators (ej)j=i_...^s. Every involution a' of Tn with a fixed point is 
conjugated with an involutive automorphism a : Zn ^ Zn stabilizing ajf : a^f o « = ajf . 
Moreover a{ei) = ±eo-(i), where a G 5*^ is a product of independent transpositions such 
that rii = na(i) Vi = 1, . . . , s. The induced action of a on 3^{Tn, K) commutes with 
A = A+ : 

(A/)oa = A(/oa) V/ G 3^(T^, iT) . 

Hence /\^{6e)\F = Ap{6e\F), where F = Fix(a). In particular, if Z^ is not harmonic 
then so is F. Choosing a appropriately, we arrive at the same conclusion as in I2.22l d. 
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2.6. Uniqueness sets. Let K he a field of characteristic 2. A subset U C vert(r) is 
called a uniqueness set for Harm^(r) if every function / G Harm^(r) that vanishes 
on U vanishes identically. Thus every (anti)harmonic function / on F is uniquely 
determined by its restriction f\U. 

The boundary of a bounded plane domain is a uniqueness set for the classical har- 
monic functions. In our discrete setting, it may happen that just a part of the boundary 
(or of the interior) serves as a uniqueness set for binary harmonic functions. Let us 
give several examples. 

Examples 2.29. 1. An extremal vertex of the linear string P„ is a uniqueness set 
for Harm^(P„). Every pair of neighborhooding vertices of the circular graph C„ is a 
uniqueness set for Harm^(C„). 

2. More generally, T x {vi} and T x {t>j,t>j+i}, 2 < i < n — 2, are uniqueness sets 
for the (anti) harmonic functions on F x P„, and F x {vi,Vi+i} is that on F x C„. 

3. The set of all extremal vertices of a finite forest F is a uniqueness set for the 
(anti) harmonic functions on F. On the other hand, the reduction T^^^ (Xred^ respec- 
tively) as in 11.21 regarded as a subgraph of F is a uniqueness set for harmonic (respec- 
tively, antiharmonic) functions on F, see the proof of 11.31 

4. Every side of a triangle n„ inscribed in a triangular plane lattice is a uniqueness 
set for IIarm^(n„). 

5. The exterior circle is a uniqueness set for the conic lattice Cn{Tn) made of m 
concentric plane copies of C„ joint one to another by radial edges, the last copy being 
also joint with a new vertex at their common center. 

2.7. Periodic harmonic extension. The idea behind 12.221 a and 12.181 is as follows, 
cf. fSuSlinEI- 

2.30. To any function / G 5'(FxC„, K) one associates a sequence /i, . . . , /« G 5'(F, K), 
where fi = /|(F x {vi}). Letting A='= = Ap^,^^ we obtain 



(A=*=/). = /(. — 1) mod n 

+ A^ (/,)+/(,+ 

1) mod ni 



,n. 



Therefore / G Harm (F x C„) if and only if 

(6) /(j+l) mod n = mod n + ^rifi) — 1' 

Starting with an arbitrary pair mq = (/o, /i) &V := [5'(F, k)]^ and applying successively 
the automorphism Jr = = ^ ^ ^± ^ G Aut{V) we extend Uq to a function / on 
F X C„ so that 

(/l, /2) =Ui = Jr(Mo) = -''r(/o, /l), • • • , (/n, fn+l) = Un = Jr(Mn-l) = ^r(/n-l, fn) ■ 

This extension / is (anti)harmonic provided that it is periodic. The latter holds if and 
only if J^{uo) = uq. By recursion we obtain 



F„_i(A±) F„(A±) 



- 1 F„(A±) F„+i(A± 
Thus 

Jr(Mo) = Mo 



F„_i(A±)/o + F„(A±)/i = /o 
F„(A±)/o + P„+i(A±)/i = /i. 

In particular (0, /i) G ker(id+jp) /i G ker(F„(A^)) n ker(id +F„_i(A^)). 
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Hence F x C„ is (anti) harmonic if and only if 1 G spec(Jp), or equivalently, if there 
exists A G spec( Jr) fl /i„. We have 

spec(J^) = {Xek"" : X + X-^ e a^} . 

Thus r X C„ is (anti) harmonic if and only if there is A G /Xn such that A + A~^ G cTp . 
This proves the second assertion in I2.22[ a. 

Remark 2.31. If F x C„ is harmonic then, according to the symmetrization and unique- 
ness principles, every (anti)harmonic function on F x C„ with /o = /i is necessarily 
symmetric i.e., fi = /(-i+i) mod n for alH = 1, . . . , n. 

We let m = card(vert(F)). Since ^ spec(Jp ), G GLfc(2m) has finite order i.e., 
3n G N : J" = 1 F„_i(A^) + id = = F„(A^). The previous discussion 

leads to the following result. 

Proposition 2.32. Vn > 3, d^{T x C„) < dimfc(K) = 2 card(vert(F)). The equality 
holds if and only if n = mod ord . 

Examples 2.33. 1. For F = C3, ord(A+) = 2 and ord(J+) = 6. So (^^(Tg,^) < 6 
Vn > 1 and (i"'"(T3 „) = 6 <^=^ n = mod 6. The cyclic shifts in the vertical 
direction of the harmonic patterns 

/l 1 1 l\ /O 1 1 1 l\ 

;ii= 1 1 1, /i2= 1 1 1 

\o 1 1 ly \o 1 1 1 oy 

form a basis of Harm'^(T3,6). The 2-sheeted covering vr : ¥35 T33 yields 4- 
dimensional subspace 7r*(Harm'''(T3 3)) C Harm'''(T3 g). 

2. Likewise, for F = C5, ord(A+j = 3 and ord(J+) = 15, so d+iT^^^) < 10 Vn > 1 
and d'^iT^^n) = 10 <^=^ n = mod 15. The 3-sheeted covering vr : T5 15 — * T5 5 
gives rise to 8-dimensional subspace 7r*(Harm"''(T5^5)) C Harm"*" (¥5^15). 

2.34. Similarly, for the path P„_i, the function / = (/i, . . . , fn-i) G J'lF x Pn-i, K) 
is (anti) harmonic if and only if © holds for all i = 1, . . . , n — 1 with fo = fn = 0. 
By recursion, the latter holds if and only if fk = Ffc(Ap)/i Vfc = 2, . . . , n. Thus 
fl G 3'{T,K) extends to / G Harm±(F x Pn-i,K) if and only if F„+i(A^)(/i) = 0. 
Hence Harm^(F x Pn-i,K) = ker(E„(A^)). This shows ITTTKl b. 

2.8. Doubling the periods. 

2.35. In this subsection we consider an injective /c-endomorphism 6 : 3'(Z^,A;) — >■ 
9^(Z^, k), which sends a function 



r 


s 


t ■■■ 


u 


V 


w 


X 


y 


z 
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into 



/ 



r 


r + s 


s 


s + t 


t 


r + u 





S + V 





t + w 


u 


U + V 


V 


V + w 


w 


U + X 





v + y 





w + z 


X 


X + y 


y 


y + z 


z 



\ 

For instance, 6 sends the harmonic function 



... 1 


1 


1 


1 


1 ■■■ 


... 1 


1 


1 


1 


1 ■■■ 


... 











■■■ 


... 1 


1 


1 


1 


1 ■■■ 


... 1 


1 


1 


1 


1 ■■■ 



to the harmonic function 



( 



... 1 





1 





1 ■■■ 


... 1 





1 





1 ■ •■ 


... 











■■■ 


... 1 





1 





1 ■■■ 


. . . 1 





1 





1 ■ ■ ■ 



replacing the bi-period n(/io) = (1,3) by the bi-period n{b{h^^ = (2,3). 

Proposition 2.36. The endomorphism 6 stabilizes the subspace Harm'''(Z^) and gener- 
ically doubles bi-periods. More precisely, if f & Harm'''(Z^) is a bi-periodic binary har- 
monic function with bi-period n{f) = {ni,n2) different from a shift of or ^Hq, then 
77-1, n2 > 1 and 6{f) G Harm'''(Z^) has bi-period n{6{f)) = 2n{f) = {2ni,2n2). 

Proof. It is easily seen that 6{f) is harmonic if so is /. If / 7^ is constant in vertical or 
horizontal direction then it is a shift of one of the functions Hq or */io, which has been 
excluded. Thus ni,n2 > 1, the function / on is non-constant in vertical direction, 
and so 6{f) possesses a nonzero line 

(■■■ u -\- X V -\- y w -\- z ■■■). 

If (mi, 0) is a period of 6{f) then necessarily mi is even and (mi/2, 0) is a period of /, 
so (mi,0) = (2^1,0) is a minimal such period. By symmetry n{6{f)) = (2^1,277,2), as 
stated. □ 

Remarks 2.37. 1. Actually 5 provides linear injections Harm'''(T„j^„2) Harm'''(T2ri,^,2n2y 

2. Doubling of just one of the periods ni,n2 is impossible in general, as Example 
12.331 1 above shows. 



^Respectively, */io to (5(*/io) = *^(^o) 
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3. A similar doubling is equally applied in higher dimensions, over any field K of 
characteristic 2. Namely, for any / G 5'(Z*,i^) and m G Z** we let 6{f){u) = f{v) if 
u = 2v has all coordinates even, otherwise 

2tiSneighb(ii) 

where 2v G neighb(-u) iff f G and the coordinates of m — 2f are equal to or ±1. 

2.9. 'Lights Out' game on graphs. The game 'Lights Out' on a finite graph F 
consists in the following Pej, (Su4j . Each vertex of F can be in one of the two states 
'on' or 'off'. A move consists in changing the state of a vertex and, simultaneously, 
of all its neighbors. The goal is to get, after a sequence of moves, all states 'off'. An 
initial position will be called pattern. A pattern is winning if there exists a sequence of 
moves terminating at a zero pattern. The graph F is called winning if the game on F 
wins starting with an arbitrary pattern. Most of the following results are well known, 
see e.g., |Su5j . 

Proposition 2.38. (a) A finite graph F is winning if and only ifV is not harmonic. 

(b) Every nucleus ofT yields a linear relation for the functions {a^}^evert(r) 

veN+ 

where 5^ is the characteristic function of N and a+ is as in (0). 

(c) The space of winning patterns Vr = span (a+ : v G vert(F)) in 3^{T,k) is the 
orthogonal complement to Harm'''(F) w.r.t. the standard bilinear form {x,y) on 
3^(F,fc). 

(d) Every antiharmonic pattern on F is winning. 

(e) For every T, the all-on pattern is winning. 

Proof, (a) A pattern can be considered as a binary function on F. The move at a 
vertex v corresponds to the shift by in 5'(F, k). Thus the game on F is winning if 
and only if the group of translations generated by (t^+ : v G vert(F)) acts transitively 
on 3'(F, k). The latter holds if and only if the functions a^, where v G vert(F), span 
5'(F, k), if and only if the matrix / + adj(F) of Ap with columns (a^) is non-degenerate, 
or, equivalently, ^ spec"'"(F), as stated. 

By definition, (5jv £ Harm'''(F) if and only if is a nucleus of F. This yields (b). Now 
(c) follows from (b). The proper subspaces Harm"'"(F) and Harm~(F) of the laplacian 
Ap being orthogonal, (d) follows by virtue of (c). Further by (c), 1 G Vr if and only 
if 1± Harm"*"(F). The latter holds indeed because for any h G Harm'^(F), the nucleus 
h~^{l) being an odd graph, by the handshaking theorem it has an even number of 
vertices. □ 

Remarks 2.39. 1. In view of l2.38l c the harmonic functions on F are linear invariants 
of the game 'Lights Out'. That is, 

{h, f + at) = {h, f) yv G vert(F), V/ G 9^(F, k), Vh G Harm+(F) . 

By virtue of 12.381 a. F is winning if and only if it does not admit a nonzero linear 
invariant. 

2. I2.38[ e is Sutner's Garden-of-Eden theorem |Su5j . The desired transformation of 
the all-one pattern into the all-zero one is achieved via moves at the vertices of an 
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odd- domination subgraph N of T. The latter means that every vertex f of F must have 
in an odd number of neighbors including v itself ii v E N, so that 1 = J2veN '^t ^ 
Vr- Given any graph F, Sutner's theorem actually proves the existence of an odd- 
domination subgraph of F (see also |Anj ). 

3. Counting harmonic toric 2-lattices and counting points on an 

elliptic cubic curve 

3.1. Constructing harmonic tori from polynomials. 

Proposition 3.1. Given a polynomial p{x) G k[x] with p{0) = 1 and a root z E¥g of 
p, one can construct a harmonic toric lattice Tn, and every such lattice appears that 
way. 

Proof. Indeed, z can be written in a unique way as z = C + C"^? where C, G Fq2. There 
is a unique decomposition 

p = 1 + 

0<ai<a2<...<as 

and p{z) = yields 

0<ai<ci2<---<"s 0<cJi<a2<---<"3 

Thus letting Xi = ("^ gives a solution of (U)) with 

n 

Hi = Old Xi = — , i = 1,. . . ,s, 

gcd(n, ai) 

where n = ordC = iordz. 

In this way we obtain all harmonic toric lattices T^^. Indeed, given n = {ni, . . . , rig) 
and m = lcm(ni, . . . ,ns), we let q = 2-^*^™'\ and we fix a primitive {q — l)-st root of 
unity ( G /ig-i. Any solution x = (xi, . . . , Xg) G {¥gY of ([T]) can be written as Xi = 
where a^nj = mod (g — 1), i = 1, . . . , s. Letting z = ( + G ¥g, the first equation 
in dH) is equivalent to p{z) = 0, where p = 1 + J2i=i ^ ^[^]- '-' 

Remark 3.2. Letting above ( = C,'^, where gcd(c, g — 1) = 1, we obtain z = Tc{z') and 
Xi = where = ^ + G Fg is a root of the polynomial Pc = I + J2i=i Tcai- 

Example 3.3. For every n = 2s + 1 odd and for every s-tuple n = (n, . . . ,n), the 
hypercubic toric lattice is harmonic. Indeed, if G A; is a primitive n-th root of 
unity then xj = , j = 1, . . . , s, gives a solution of (HJ with nj = n^j. In particular 
T(5^5), T(7_7 7), T(ii 1111 1111) etc. are harmonic. Therefore bv I2.22l d. so is Tn provided 
that Hi = mod 5 for at least 2 values of i, or = mod 7 for at least 3 values of i, 
etc. However, C5, T(7 7) and T(ii_ii) are not harmonic, see 15.11 a. 13.71 and Appendix 1 
below. 

3.2. Partners. 

3.4. Let E be the affine plane cubic with equation 
(7) {l + x + y){l + xy) = l, 

and set E* = E \ {(0,0)}. By virtue of 12.121 b the toric 2-lattice T^.n is harmonic if 
and only if the curve E*{k) possesses a point {x,y) with ordx | m and oidy \ n. 
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We are interested in the infinite table £ composed of all pairs (m, n) G such that 
the lattice Tm,n is harmonic, or equivalently, such that x™" = = 1 for some point 
{x,y) G E*(k). We consider also the subtable 

So = { (ord X, ord y) : {x,y) E E*{k)} C 8, . 

Thus £o is the set of all bi-torsions of points on E*, or in other words, the set of all 
minimal bi-periods of double periodic binary harmonic functions on 7?. Notice that £o 
contains a set of primitive generators of £ viewed as a module over the multiplicative 
semigroup (cf. 12. 221 d). We call m and n partners if {m,n) G £o. For instance 
(47, 178481) is a pair of partners found by Zagier. 

The following lemma is a reformulation of Theorem 5.2 in |HMPj : the latter also 
covers the case of square grids over F3. 

Proposition 3.5. The table of partners £0 is infinite. 

Proof. If £0 were finite there would exist a prime p such that 

p > M = max{m : (m, n) G £0 for some n G N} > 5 . 

We have d = 1 mod p Vrf | (2^ ^1)- Indeed, for every prime divisor Z of 2^ — 1, 

2^ = 1 mod / =^ ord; 2 = p =^ I = 1 mod p . 

It follows that d = 1 mod p and so, (i > p > M if (i > 1. 

However the toric lattice Tg_i being harmonic bv 12.131 we must have {q — l,q — 
1) = {kdi,ld2) for some di, d2, k, I E N such that ((^1,^2) ^ £0 di > 1- Since 
di > M this yields a contradiction. □ 

Somewhat more precise information can be deduced by using the Hasse-Weil formula. 

3.3. Hasse-Weil formula. The cubic curve E as in 13.41 has 3 points at infinity: (1 : 
: 0), (0 : 1 : 0) and (1:1: 0). Hence the projective closure ^ of is a smooth 
elliptic curve. For g = 2*" we let Sr {sr = 3^ — 4, respectively) be the number of points 
on E{¥g) {E*{¥g), respectively). 

Lemma 3.6. We have Sr = q {l — (a+ + «1)) — 3, where a± = (—1 ± \/—7)/A are the 
complex roots of the polynomial 2t^ + t + 1. Moreover, the Hasse inequalities hold: 

(v^-l)'<s, <(v/^+l)^ 

Proof. Since si = 4, the Hasse-Weil formula |Ko[ Ch. V, §1, Exercise 7] gives in our 
case: 



00 



j2 'ff = log <M = log j^^^0^^ = 2: ^ (1 + 2'(i - («; + . 

Now the assertions follow easily. □ 

From ESI and |HMP| §5, Remark] we deduce the following results. 

Corollary 3.7. (a) Vg = 2'' > 16, ^*(F,) ^ 0. 

(b) All Mersenne primes g — 1 = 2^ — 1 with p > 3, and all Fermat primes q+1 = 
2^+1 with I > 1, are self-partners i.e., (g±l, g±l) G £o- Whereas for p = 2,3 
one has (3, 3), (7, 7) £9. 
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Proof, (a) follows by virtue of 13.61 

(b) According to l2.13t the toric lattice Tg±i^g±i is harmonic Vg = 2^, except for Ti i 
and Tjj. Thus for all those m = n = q ± 1, Q has a solution If g ± 1 is prime 

and ^, ?7 7^ 1 then ord^ = ord?7 = q ± 1 that is, (g ± 1, g ± 1) G £0- This proves the 
first assertion. 

Since x = 1 or y = 1 for every point {x,y) G E*{A) then (ordx,ord?/) = (1,3) or 
(3, 1) and so (3, 3) Sq. Neither (7, 7) G £0 since S3 = 0. □ 

So far only 5 Fermat primes and at most 43 Mersenne primes were found, see e.g., 
jLLMPl 7.3], [Mil, M- 

The computer findings in Appendix 1 suggest the following conjecture, cf. 12.131 

Conjecture 3.8. Vg = 2^ (r > 6), g — 1 and g + 1 are partners and auto-partners that 
is, (g ± 1, g ± 1) G £0 and (g + 1, g — 1) G £0. 

The latter does not hold for r = 5. Indeed (31, 33) G £ \ £0, see 12.131 and Appendix 1. 

Examples 3.9. 1. (5, n) G £ if and only if n = mod 3 or n = mod 5. In particular 
(5, 5) G £0 is a self-partner, and there is no further partner of 5. 

Indeed, as (1,3) G £ then {k,3l) G £ for every k,l > 1, in particular for k = 5. 
Further, 5 is a primitive self-partner since for any primitive 5-th root of unity ( E fi^, 
the pair (x, y) = {(, satisfies ((H) with s = 2, m = n = 5. Consequently by virtue of 
EMs^, {5k, 51) G £ VA;, / > 1 and so (5, 5/) G £ V/ > 1. 

Conversely, if (x, y) G E* and (ordx, ordy) = (5, n) then x G Fig \F4 satisfies = 1 
and y = x"^, satisfy 

y"^ + {z + l)y -1-1 = 0, where z = x + x~^ . 

Thus n = oidy = 5. 

2. (7, ra) G £ if and only if n = mod 3, and 9 is the only partner of 7. 

Indeed, (7, 3A;) G £ VA; > 1 because (1,3) G £. If {x,y) G E* and (ordx,ord?/) = 
(7,n) then x G Fg \ F2. We have (7,n) G £0 =^ fo{n) = /o(7) = 3 =^ n\ 
(2^ ± 1) =^ ne {7, 9}. But n 7^ 7 as S3 = and so, E*{¥27) = 0. Hence n = 9. 

Since (1,3), (7,9) G £0, the latter set properly contains the set of all primitive 
generators of £ over N^. 

3.4. Partnership graph. We observe that: 

• For every (m, n) G £0, both m and n are odd (cf. 12.221 c). 

• Every odd n G N has a partner, and the number of these partners is finite. 

• (m, n) G £0 <^==^ {n,m) G £0. 

Thus the partnership defines an equivalence relation on Nodd- Answering a question of 
the author, Zagier proposed the following 13.1 01 ITTT] and ITT^ below. Our proof of 13.111 
based on 15.91 is somewhat different from the original one. 

We let below div(n) (div*(n), respectively) be the set of all (proper) divisors of 
n G N. For g = 2'' we write for short div(g ± 1) meaning div(g — 1) U div(g -|- 1). 

Definition 3.10. We let T^^^ be the infinite graph with loops such that vert(T*^^)) = 
Nodd and [m,n] G edg(T^^)) <^==^ {m,n) G £0. We call 7^^^ the partnership graph. 

Theorem 3.11. All connected component of?^^^ are finite. 
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Proof. We let Vr be the subgraph of J'''^^ with vertices in the finite set 

Vr = {ne Nodd : n \ (2'^±1)}. 

Given n E Nodd, we let T^^^n) be the connected component of 7^^^ which contains the 
vertex n. We claim that the function fo{n) is constant on each connected component 
of y^^^. In particular T^^^n) C Vj^^^n)- The level sets Vr of /o being finite, this proves 
the theorem. 

To show the claim we note that, due to 12.121 a. [m, n] G edg(y*-^'') if and only if 
^ + ^^^ = 1 + ri + ri^^ for some primitive roots ^ G fim and rj G /in- According to lS.lUl a. 

/o(m) = deg(C + C^) = degir] + r]~^) = fo{n) 

and so, the claim follows. □ 

Notation 3.12. We denote by S{m, n) the set of all solutions (^, rf) of ((Tj) of type (m, n) 
that is, with ^ G Hm (respectively, rj G /in) being a primitive m-th (respectively, n-th) 
root of unity. We label the edges [m, ra] G edg(T*^^^) with s{m,n) = | card(S'm,n) G N, 
with one exception: instead of the edge [1, 3] we introduce two directed edges, [1 3] 
labeled by 1 and [3 1] labeled by 2. Clearly, 

Sr = 2 s(m,n), VrGN. 

m,nGdiv{g— 1) 

Moreover this labeling possesses the following properties. 
Proposition 3.13. \/n G Nodd and Vg = 2^ r > 3, 

(a) EmeN^dd ='^(^) ^■ 

(b) EdGdiv(n),mGNodd^(^'"^) 

(d) 2 E... 

(e) s(g - 1, g - 1) + s(g + 1, g + 1) + 2s{q - 1, g + 1) > 2 (^(g - 1) + (^(g + 1) - g). 

Proof, (a) holds because for every {CiV) ^ S{m,n), the pairs {rj,rj~^) and (C, C~^) 
uniquely correspond to each other. Since Ed|n '^('^) ~ ^) follows from (a). Summing 
up (a) over the edges of yields (c). It is easily seen that d G div(g ± 1) <^=^ 
/o(c?) I r <^=^ d eVs for some s G div(r). Hence 

^^¥;(n)= ^ ^(n)= ^ ^{n)+ ^ </^(n') = 2g. 

s|r ngVs nediv(<j±l) n|(<5r— 1) n'\(q+l) 

Thus the summation of (c) over the set div(r) yields (d). 

(e) By virtue of 13. 131 a. s{d, g — 1) + s{d, g + 1) < '{>{d). Moreover 

¥?(g - 1) + ¥?(g + 1) = s(g - 1, g - 1) + s(g + 1, g + 1) + 2s(g - 1, g + 1) 
+ (s(rf,g-l) + s(c?,g+l)) . 

(iGdiv*(g±l) 

Hence 

s(g - 1, g - 1) + s(g + 1, g + 1) + 2s(g - 1, g + 1) > L(g - 1) - ^ ^{d) 

\ dediv*(q-l) 



Hereafter if stands for the Euler function. 
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+ V^(g + 1)- Yl V{d') \ =2i^{q-l) + ^{q + l)-q) . 



\ d'gdiv*(g+l) 



□ 



Remarks 3.14. 1. If {Crj) is a solution of (|7j) of type {171,11) then {C'^i'rf ) is as well 
such a solution. Thus the Galois group Gal(F(^)) acts freely on S{m,n) and so, its 
order f{n) divides s{m,n). Letting k{m,n) = s{m,n)/ f{n), from l3.l51 a we obtain the 
equality 



Indeed, there are v?(n) = f{n)g{n) primitive n-th roots of unity. 

2. For q = the inequality ip{q — 1) + ip{q + 1) > g does not hold in general, 
although it holds at least for all r < 150 (a MAPLE checking)^. Moreover, according 
to A. SchinzeF, ^9(2^' — 1) + V5(2'' + 1) can be << 2'^/(loglogrt)^/^'^ for infinitely many 
r, although ipi^!^ — 1) + ip{2'' + 1) > 2*^/ log log r Vr G N. This can be seen as follows. 
Let Tt be the least common multiple of all numbers (p — l)/2, where p runs through 
primes = 3 mod 4 less than t. By Euler's Theorem and the quadratic reciprocity law, 
p divides 2*^^"^^/^ — (— hence 2*"* — 1 is divisible by the product of all primes 
= 7 mod 8 less than t, while 2''* + 1 is divisible by the product of all primes = 3 
mod 8 less than t. Using the Mertens formula for primes in arithmetic progressions, 
we obtain that 93(2^* ± 1) << 1/logt << l/(loglogrj)^/'*. A slight modification of 
the argument increases the exponent 1/4 to 7/24, which gives the claim. But a gap 
between the lower and the upper bound remains. Inded from a theorem of Erdos |Erj , 
and from the inequality between phi and sigma functions (HW^ Theorem 329] it follows 
that ¥?(2^ - 1) + (/?(2^ + 1) »271oglogr. 

3. Similarly to one might consider an infinite hypergraph T = Us>iT*^*-' with set 
of vertices Nodd such that (no, . . . , iis) is an s-simplex of T*^**^ if and only if (no, . . . , n^) 
is the multi-order of a point on the affine hypersurface as in (with s replaced by 
s + 1). Evidently, every s-tuple (ni, . . . , n^) G N^^^ is a face of an s-simplex in J**^*). 

4. Appendix 1: Connected components of the partnership graph 

The first 13 connected components of the partnership graph J"*^^) are shown below. 
They were found by Zagier with PARI. The labeling of the edges is according to 13.121 
We recall (see the proof of I3.1H) that the value of /o(n) equals r for every vertex n 
of V^, and this determines V^. The value of /(n) equals r if n is not underlined, and 
2r otherwise. The edges [m,n] (the loops [n, n], respectively) correspond to harmonic 
toric lattices „ (T„_„, respectively) with (m, n) G £0. 



Computations with Pari/GP done by Gottfried Barthel (a letter to the author) confirm the in- 
equality in the range r < 275. 
letter to the author. 




m 
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We observe that the graph V12 is not planar, in contrast to V,. with r < 11. 
These computations suggest the following 

Conjecture 4.1. Vj. is connected Vr 7^ 5. In other words, the connected components 
of 7^^^ are V^. for r 7^ 5 and the two components of V5. 

5. Appendix 2: Chebyshev-Dickson and Fibonacci polynomials 

5.1. Chebyshev-Dickson and Fibonacci polynomials. These polynomials T„, En 
and Fn provide an important tool for analysis of harmonicity. Indeed, as we have 
seen in 12. 2| T„, respectively. En is the characteristic polynomial of the laplacian A^^, 
respectively, Ap^, where C„ stands for the circular graph with n vertices and P„, denotes 
the path of length n. We give an account of some of their properties in 15. Ill5^ below 
according to |LMT[ Ch. 2], |Su8j . |BR,j . ICjKWj . see also references therein. 

Definition 5.1. Consider the linear recurrence 

(8) pn+i = xpn +Pn-i, where pi G k[x] Vz > . 
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Thus 

Pn ^ ( PO 

Pn+1 J \1 X J \pi 

The Chebyshev-Dickson polynomials of the first, lespectivelj, second kind T„, En G k[x] 
and the Fibonacci polynomials Fn € k[x] are defined via (jHI) by the initial conditions 

To\ /0\ f f l\ f Fo\ f 



, - . ,, respectively. 
Ji/ V^/ V-^i/ V^/ V-^i/ V-'^ 

Thus degT„ = deg-En = degF„+i = n, the polynomials T2n, E2n, -F2n+i are even and 
T2„+i, -E2n+i, -p2n are odd. They are related via 

(9) Tn = xFn = xEn-i, whcrc E^i = . 

So any property of one of the sequences (T„), (En), (Fn) is enjoyed by the other two 
up to evident changes. Notice that Fn{l) is the n-th Fibonacci number modulo 2. 

The following identities hold, see e.g. [RH ILMT] . |BR1 §4], |WP] . 

Proposition 5.2. Vm, n G N, Vg = 2*", r >1, we have 

(a) = z{z'^ + xz + is the generating function of the sequence (Fn). 
(b) 

i=0 V ^ ^ j=0,...,n, jEEn mod 2 

where ?7„(cosx) = stands for the Chebyshev polynomial of the second kind 
overR. ^° 

(c) Fg^i + -Fg+i = xFq = x'^ . Furthermore, 

=x'^ + = I 1 + 

(d) G 

= 2Tr2^(;2-i), Fg+i{z) = z{l + Tt^^{z-^)Y and Fq_i{z) + Fg+i{z) = z. 

(e) Em+n = EmEn + Em-lEn-1- 

(f) E2n = xEnEn-i + l = El + El_„ E2n+i = xEl and Zti ^2. = El 

(g) Vn > t > 0, Tn+t + Tn-^t = TnTf. In particular for n = t mod 2, + Tf = 

F ^n-{-t F ^n — t , 
2 2 

(h) FmOFn = FnOFm = F„ 



r-1 
j=0 



- mn ■ 



(i) T,„(x)=T„(x'^) = r„«(x). 



Proof. If r = P„ is a linear graph with n vertices and e = n — 1 edges then (|)^ = {"'~^) ■ 
So (b) follows by virtue of II. 11 a and 12. 2( whereas ll.ll c yields (e). In turn (e) implies 
(f), (g). Now (c) follows by recursion and implies (d). The assertions (a), (h) and (i) 
can be deduced, by virtue of (b), from the analogous identities for the usual Chebyshev 
polynomials over M, see e.g., |LMT| IKi] . □ 



"'^'^ Attention: our enumeration of classical polynomials does not coincide with those used in MAPLE. 
It is so chosen in order to write the identities in a more elegant way. 
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5.3. As before, ord^ G Nodd denotes the multiplicative order of an element ^ & . 
In the next proposition we indicate certain divisibility properties and factorization of 
the Chebyshev-Dickson and Fibonacci polynomials according to |LMT| Ch. 2], |Su3j . 
[HKWj . ^BR. §4] and ,WPJ . 

Proposition 5.4. (a) Vn > 0, G A;^, T„(^ + ^-^) = ^" + 

(b) G A;^, ord^ = min{n > : T„(^ + ^-^) = 0}. 

(c) Vn G Nodd o,''T'd for any primitive n-th root of unity C G we have 

(n-l)/2 

T„(x)=x n {x + C + (-'?■ 

i=l 

Consequently, z = ^ + runs over the roots ofTn when ^ runs over fin. 

(d) E„_i(0) = n = Omod 2, = n = Omod Banc? 
(x^ + X + 1) I En^i <^=^ n = mod 5. 

(e) Vm,n G N, gcd(Tm,T„) = rgcd{m,n) and gcd{Em-l, En-l) = -Egcd(m,n)-1- 

(f) Td I T„ I <^=^ d I n. 
5.2. Irreducible factors of Fibonacci polynomials. 

5.5. Every z ^ k can be written in a unique way as z = ( + where ( and are 
the roots of fz{x) = + zx + l G k[x\. By virtue of l5.4[ c. every irreducible polynomial 
r G divides one of the F„ |Su3| 3.1]. Namely, if r 7^ x and + ^~^) = then 
r I Fn with n = ord^. 

Remark 5.6. An element ^ G {q = 2^) can be written as 2; = (m^ + -u)"^ for some 
M G Fq \ F2 if and only if it can be written as z = C, + for some ^ G F^ , see I5.21 d 
and I5.4[ c. 

These observations lead to the following definition |GKWl ISu3j . 

Definition 5.7. Letting Irr[a;] be the set of all irreducible polynomials in k[x], we 
remind that for r G Irr[x], ordr = ordz G Nodd, whenever 2; is a root of r. For t ^ x 
we define its Fibonacci order 

fordr = ford 2; = min{n > : r | F„} = ordC G Nodd, where 2; = C + C ""^ • 

5.8. For n G Nodd, — 1 = Hdin ^ product of cyclotomic polynomials 

(10) $d(x)= n ^= n (^-^^)' 

T6lrr[a;], ordr=d l<i<d— 1, gcd(i,ii)=l 

where C G /i^ is a primitive (i-th root of unity. Hence deg$rf = v^(n), and $d is a 
product of gid) = ip{d)/ f{d) distinct irreducible factors of the same degree f{d) |LNt 
2.47]. By virtue of (fTUj). is auto-reciprocal, that is $^ = where 

* : p(x) I — > x'^''^'^p{x^^) 

is an involutive automorphism of the multiplicative semigroup 

U[x] = {pe k[x] : p{l) = 1} . 

""^^By making use of l5.2l h. one can deduce in the same way the roots of r„ for any even n. 
^^fordr is called Fibonacci index of r in |GK W| . and depth of r in |Su3| . 
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It occurs that either all irreducible factors of are auto-reciprocal or none of them 
is, depending on d. More precisely, the following happens. 

Lemma 5.9. Let ( ^ fi^ be a primitive root of unity of odd order d. We denote by 
'^i^) — 11^=0 C^^) ^ k[x] its minimal polynomial of degree r = f{d). The following 
conditions are equivalent: 

(i) r is auto-reciprocal, or palindrome. 

(ii) 2^o('^) = -1 mod d. 

(iii) degC = 2deg(C + C-'). 

(iv) Tr,(^^.^-i)(C + C-^)-i = l. 

Proof. T being irreducible of degree > 2, we have r(l) = 1. But 1 is the only fixed 
point of the involution z i — > on . In case that t = t* this involution acts on 
the roots of r, hence f{d) = degr is even. 

The roots of r being the conjugates of ( in k, t is auto-reciprocal if and only if ( 
and are conjugated. By virtue of 12.11 this yields the equivalence (i)<^==^(ii). The 
condition (iii) holds if and only if the polynomial fz as in 15.51 above is irreducible, 
where z = ( + Thus the equivalence (iii)<^==^(iv) follows, see e.g. |McEt 8.13]. To 
show the remaining equivalence (i)<^==^(iii) we consider the Laurent polynomial g{x) = 
r(x)r(x^^) G k[x,x^^]. It is auto-reciprocal that is, g{x) = g{x^^) or, equivalently, 
g{x) = h{x -\- x^^) for some h G k[z] of degree r. As h{( -\- (^^) = g{() = 0, (iii) holds 
if and only if h is reducible. 

Supposing (i) we have r = 2s, where s G N, and r(x^^) = x~^t{x). Thus h{x-\-x^^) = 
{x~^t{x))'^ =: g'^{x). The Laurent polynomial g being auto-reciprocal, it follows that 
g{x) = h{x -\- x^^), where h G k[x], degh = s and h = h^. Clearly, h is the minimal 
polynomial of C + C^^- This yields (iii). 

Conversely, let (iii) holds i.e., h = /11/12 is reducible, where hi G k[x], deg/ij = Tj > 1, 
i = 1,2 and ri -\- r2 = r. Letting gi{x) = x^^hi{x -\- x~^) G k[x] we have deggi = 2rj, 
i = 1,2. Furthermore 

t(x)t*(x) = x'^t{x)t{x^^) = x'^h{x -\- x^-^) = gi(x)g2{x) . 

Since r, r* G Irr[a;], up to interchanging gi and g2 we obtain gi = r, g2 = t* . Hence 
r = 2ri = 2r2 is even and 

r(x) = gi{x) = x^^'^hi{x -\- x~^) = x'^t{x^^) = t*{x) . 

Thus r = r*, so (i) holds. □ 

Corollary 5.10. (a) For every primitive d-th root of unity ( & fid we have deg C = 
f{d) anddeg{C + C-^) = fo{d). 
(b) Let t{x) = X]i=o^«-^* minimal polynomial of (, and let t{x)t{x~^) = 

1 -\- X]j=i ^ji^'' + x~^). Then the minimal polynomial ri{x) of ( + (^^^ is 

ri + E;=i^s-,T,(a;) 
ri{x) = < with degr] = 

ii+E;=i5.T,(x) 

Following |GKWj . |Su3t 3.2], |SBj we list below some important features of irreducible 
factors of the Fibonacci polynomials. 

I30r, equivalently, [¥2(0 : F2(C + C"')] = 2. 



s = r/2 if T = T*, 
otherwise . 
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Proposition 5.11. For every r G Irr*[x] with fordr = d, the following hold. 

(a) degr = /o((i) i.e., the splitting field of t is ¥q with q = 2-f°^'^\ 
(h) t\F, ^ d\n ^ Fd\F^. 

(c) d \ {q — 1) (and so r \ Fg_i) if and only if the linear term of t vanishes, if and 
only if the polynomial fz{x) = x'^ + zx + 1 E ^^[x] splits over ¥q, where z E¥q 
is a root of r. Otherwise d \ {q + 1) (and so r | -Fq+ij. 

(d) Vn = 2A; + 1 G Nodd, 

Fn= n ^' = ^L 

r Girr* [x] , ford t \ n 

where Rk = Fk+i + Fk, deg Rk = k, Rk is square-free and contains the monomial 
x''~^. The splitting field of Fn is ¥g, where q = 2-^"*^"^. 

(e) Vg = 2^ , r > 1, the splitting field of Fq^i is ¥q. Moreover 

F,-iF,^i = n = (^''' + • 

rGlrr* [x], deg T\r 

(f) For every odd prime p and for every irreducible factor t of Fp, 

ifip) if f{p)^l mod 2, 

\f{p)/2 otherwise. 



deg r 



In particular Fp = if and only if either f{p) = p — 1 or f{p) 
mod 2. This cannot happen if p = ±1 mod 8, p > 1. 



p-i 



Proof, (a) follows by virtue of l5.1(H a. and (b) follows from I5.4l b.f. If 2; is a root of r 
then by (a), 2 G F^^ and so, by virtue of Od, Fg^i{z) = Tr^^{z-^) = 

and Fq+i{z) = <^==^ TiY^iz^^) = 1. By (b), | (g — 1) in the former case and 
d \ (g + 1) in the latter one. This yields (c). As T„ = (d) and (e) follow easily 
fromlOlf andlOc. For the proof of (f), see jGKWj . □ 

Remark 5.12. We let Irri[x] = {r G Irr[x] : r = x'^''^^ + x"^^^^'^ + ...}. Thus 
r G Irri[x] if and only if r G Irr[a;] and TTk{z){z) = 1 for any root z = ( + of r. 
By virtue of (d) above, an odd number of irreducible factors of Rk belong to Irri[a;]. 
Hence for every n G Nodd, there exists d \ n and a primitive ci-th root of unity ( E fid 
such that 

/oW-i 

Tr.(.)(c + c-^)= E r + rn = i. 

1=0 

5.13. For any n G Nodd, the following analog of the cyclotomic polynomial $„ (see ITU)) 
were introduced in |Su3j : 

Pn= n 

rglrr* [x] , ford r=n 

The following properties of these polynomials were established in |Su3j . 

Proposition 5.14. (a) degp„ = ip{n). Furthermore, pn has irreducible fac- 

tors, all of the same degree fo{n), the same multiplicity 2 and with the same 
linear term. 
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(b) Vg = 2\r> 0, 

Fn = X{pd and F,n = x''-^F^ = x'^-^\[pl. 

d\n d\n 

(c) By the Mobius inversion formula, 

d\n 

5.3. +-involution. The automorphism group Aut(A;[x]) is isomorphic to Z/2Z and 
consists of the identity and the involution 

p : k[x] — ^ k[x], p{x) I — > p^{x) := p{x + 1) . 

Following [GKW| we call p'^ the conjugate of p. This notion plays an important role 
in the analysis of harmonicity of plane grids, see Section 2. The ring of invariants 
k^[x] := k[xY = ker(5), where 5 = p + id G End(A;[x]), consists of all self-conjugate 
polynomials p = p'^. Clearly, p preserves degree and irreducibility. 
The following proposition extends Lemma 15 in |(TKWj . 

Proposition 5.15. (a) im{5) = ker(5) = vect {5{x"') : n G Nodd)- 

(b) We have deg6{p) = degp — 1 G k[x] with degp G Nodd, and degp G Ncven 

G k~^ [x] . 

(c) k[x] = k^[x] © fcoddN; where fcoddN ^ k[x] is the subspace of odd polynomials. 

(d) 6{fg) = f6{g)+g6{f)+6{f)6{g), ^f,g G k[x]. 

(e) In particular S{fg) = fS{g) V/ G ker((5), Vf? G k[x]. 

(f) = (5(j9))^, Vp G k[x]. Consequently, k~^[x] = ker(5) is stable under the 
endomorphism p i — > p^ . 

Remark 5.16. 1. By (a) and 15. 221 c below, every polynomial g G k'^[x\ is of the form 

n 

for some odd p E k[x] and some gi G k[x]. 

Similarly every polynomial / satisfying /"*" = / + 1 is of the form f = x + g for some 
g G For instance, by virtue of l5.2l c this is so for the polynomials fg = Fg^i+F^-^ 

Vg = 2\ 

We let Irr"'"[x] = Irr[x] fl k^[x] C k[x] be the set of all irreducible self-conjugate 
polynomials. 

Examples 5.17. 1. The polynomials in Irr'*'[x] of degree < 8 are the following ones: 

X^ +X + 1, x'^ + X + 1, x^ + x^ + x^ + x^ + 1, 

x^ + x^ + x^ + x^ + 1, x^ + x^ + x^ + x"^ + X^ + X + 1 . 
2. F5 = {x^ + a; + 1)^ is the only self-conjugate Fibonacci polynomial |GKWj . 

5.18. For any q = 2^, we let Trq(x) = YllZo ^^^^ Tr^(x + y) = Trg(x) + TTq{y). 

We let hr{x) = x'^ + x + 1 G k[x] and hr{x) = 1 + Trg(a;). The following proposition 
can be checked readily, see |LJN| Theorem 3.80] for (b) and |LJN| Exercise 3.90] for (c). 
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Proposition 5.19. (a) Vg = 2*", r > 0, we have 

hr{x) = Kixix + 1)) = 1 + Tif^ix) + T4^(x) e k+[x] 

and 

hrix) = 1 + Trg2(x) = 1 + Tr^(x) + Tr^(x) = KiTigix)) E k+[x] . 

(b) The decompositions of hr, hr into irreducible factors overWq are, respectively, 

q/2 q/2 

hr{x) = + Pj) and hr{x) = + x + Pj) , 

i=i i=i 

where (3j G runs over the affine subspace {z E¥q : TrF,(-2) = 1}. 

(c) \/u,v G k, hr{u+v) = hr{u)+hr{v) + l. The splitting field of hr is¥g2. Moreover 
hr has simple roots that fill in the r-dimensional affine subspace 



{ze¥g2 : z'^ = z+l} = {ze ¥g2 : p{z) = pg{z)} C F„2 



-1/ 



where p : z i — >■ z + 1 and pg : z i — > z'^ . This subspace is parallel to hr (1) = ¥q 
and stable under the Frobenius automorphism and under the involutions p and 

Pq- 

Ve G N and Vm G Nodd, we let D2{2'^m) = e. 

Proposition 5.20. (a) hrohg = hgohr = hr+s+hr+hs- Consequently, h2s = hsohg. 
More generally, \/q = 2^ and Vs G N, 

hgs = hgo - yO hs^ . 

q 

(b) hs I hr <^==^ r = ms, where m G Nodd- In particular hi \ hr ^^==^ r G Nodd- 
(c) 

hgcd{r,s) if D2{r) = D2{s), 
otherwise. 



gcd(hr, hs 



Similarly, 



gcd(/tr, hs 



hgcd{r,s) if D2{r) = D2{s) 
1 otherwise. 



Proof. The proof of (a) is easy and can be omitted. To show (b) we assume that hs \ hr, 
and we let q = 2"^, q' = 2*. By virtue of 15. 191 b. ¥gi2 C ¥g2, hence s \ r. If r = ms then 
for any root z of hg, 

if m is even, 
otherwise. 



^q ^ ^q'- 



Thus hs I hr implies that m G Nodd- The converse is easy. 

(c) Assume that hs and hr have a common root z E k of degree 6. By virtue of 
15. 191 c. ¥26 C F22r n F22S, hence 6 \ 2gcd(r, s). Since z^"^ = z"^" = z + 1 z, we have 
z G (F22r \ F2r) n (F22. \ ¥2s). It follows that 5 does not divide gcd(r, s) i.e., 5 = 25', 
where 6' \ gcd(r, s). 
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Letting r = m5' and q = 2^' we obtain 



or ^71 

z + l = = 



z if m is even, 
z^ otherwise. 



Hence m e Nodd- Similarly, s = nS', where n G Nodd, and so, D2{r) = D2{s) = D2{S'). 
Moreover, hs'{z) = 0. It follows that gcd{hr, hg) \ /igcd(r,s)- Vice versa, by virtue of (b), 
hgcd{r,s) I gcd{hr, hs). Therefore (c) follows. □ 

5.21. For a polynomial t = + ar-ix^~^ + . . . G Irr[a;] we denote Tr(r) = a^-i = 
Tt^^{z), where q = and 2; G is a root of r. We let 



Irrjx] = {r G Irr[x] : Tr(r) = i}, i = 0, 1, 
and 

Irr^[x] = {tt^ : r G Irr[,T] \ Irr^[x]} . 

Proposition 5.22. (a) Every irreducible factor of hr belongs to Irr^[x], and every 
T G Irr"'"[a;] divides one of the hr, r > 1. 
(b) Every polynomial r G Irr^[x] of degree 2r admits a decomposition 

r-l 



i=Q 



where p e¥g (q = 2") and TrF,(/5) = 1. 

(c) The map 

a : k[y] k'^[x\ = k[x{x + 1)], q{y) 1 — > p{x) := q{x{x + 1)) , 

is an isomorphism of k-algebras. 

(d) Irr(A;+[x]) = a(Irr[x]) = Irr^[x] |J Irr~[x] . Moreover a {In i[x]) = Irr'''[x] and 
a(Irro[x]) = Irr^[x]. 

Proof. An element z G A; is a root of one of the polynomials hr, r > 1, if and only if z 
and z + 1 are conjugated, if and only if the minimal polynomial r of 2 is stable under 
the involution p : z 1 — > z + 1 i.e., r G Irr'*'[x]. Hence (a) follows, 
(b) is immediate from l5.1^ b and I5.19l b. 

To show (c) it is enough to establish that a is surjective. For every / = HILi '^i ^ 
k^[x] of positive degree, the involution p acts on the set {Tj}j=i .. „ of all irreducible 
factors of /. Since f = f^ = nr=i ''"i^ then either Tj G Irr^[x] or = tj 7^ Ti 
Vi = 1, . . . , n. 

For every f = 111=0 + z"^') G Irr[x], 

r—1 r—1 

{ff+){x) = l[{x + z^'){x + 1 + z^') = l[{x^ + X + 13^') = a{T){x) , 

i=0 i=0 

where {3 = z{z + 1) G Fg, g = 2*", and r G k[x\. 

We will show below that for every r G Irr^[x], t = a{f) for some f G Irr[2;]. Thus 
/ G k[x{x + 1)]. Hence k^[x] = k[x{x + 1)]. Now (c) follows. 
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(d) Let [3 G Fg, where q = 2*", be a root of a polynomial r G Irro[x] of degree r > 1. 
Then TrF,(/5) = 0, hence /? = + 1) for some z e¥q jLNl 2.80]. Therefore 

r— 1 r— 1 r— 1 

j=0 i=0 j=0 

where f{x) = YYjZoi^ + z"^^) G Irr[x] as degz = r = deg f. 

Suppose further that f = f^. Then = 2; + 1 for some zG{l,...,r — 1}. Hence 
= /?, and so deg P < i < r, a. contradiction. This proves the last equality in (d) (cf. 
[TM Exercise 3.86]). 

For every r G Irr^[x] of degree degr = 2r > 0, the involution p acts on the set 
{zj = z^^ }j=o,...,2r-i of roots of r, where z is one of these roots. If p{z) = z+1 = z^' = Zi 
then p{zj) = Zj + 1 = Z(^i+j) mod 2r, j = 0, . . . , 2r — 1. As = id we have j + 2i = j 
mod 2r Vj, hence i = r. Therefore 

2r-l r-1 

^(2^) = n + = + + -^i + 1) = t'{x{x + 1)) = a(f)(x) , 
i=o j=o 

where 

r— 1 r— 1 

f{x) = Hix + (3j) = Hix + f3'') G k[x] 

j=0 j=0 

with Pj = Zj{zj + 1) and [3 = z{z + 1) G Fg, q = T' . Since deg/5 = r = (deg2;)/2, we 
have f[x) G Irr[x]. Moreover Tr]f_j(/?) = 1 since the polynomial 2;^ + 2; + /3 G Fq[x] is 
irreducible over Fg. Hence f G Irri[a;]. This shows that a(Irri[x]) = Irr'^[x]. Now (d) 
follows. □ 

Remarks 5.23. 1. ( |LN1 Exercise 3.87]) For any odd prime p such that ordp2 = p — 1, 
'^{^) = ^^j^ ^ Irri[x] and so, a(r) G Irr"^[x]. 

2. K e Irr+[x] r G {1,2}. 

3. Vr = 2", /ir has 2^'~"~^ irreducible factors of the same degree 2r = 2""^^. 

4. We have hi{x) = 1 + Tr24(a;) = rr+, where r, r+ G Irr[x], r = + x'^ + 1 7^ 
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